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ABSTRACT 

The recently published GALFA-HI Compact Cloud Catalogue lists 20 neutral hydro- 
gen clouds that might pinpoint previously undiscovered high-latitude dwarf galaxies. 
Detection of an associated gamma-ray dark matter signal could provide a route to 
distinguish unambiguously between truly dark matter dominated systems that have 
accumulated neutral hydrogen but have not successfully ignited star formation and 
pure gaseous structures devoid of dark matter. We use 4.3 years of Fermi observations 
to derive gamma-ray flux upper limits in the 1-300 GeV energy range for the sam- 
ple. Limits on gamma rays from pair annihilation of dark matter are also presented 
depending on the yet unknown astrophysical factors. 
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1 INTRODUCTION 

The number of dark matter subhalos surrounding the Milky 
Way has puzzled dark matter aficionados for over more than 
a decade. While cosmic large-scale structures are well de- 
scribed by the Lambda Cold Dark Matter (ACDM) cos- 
mological model, a number of discrepancies exist between 
the standard theory of galaxy fo rmation and observat ions 
of substructures at smaller scales IjFrenk fc Whitell2012l ). In 
particular, numerical ACDM simulations consistently pre- 
dict that galaxies must be sur rounded by a huge popu - 
lation of subhalos (iKlypin et al. 1999; M oore et all 1 1999? ) . 
Intuitively, subhalos would encompass anything from the 
largest satellites {e.g. dw arf galaxies) to substructu res with 
masses around 10"'' M© (jLoeb fc Zaldarriagal 12003 ). Unfor- 
tunately, completing a survey of the smallest gravitationally 
bound systems is not straightforward ( Ando et al. 2008 ). At 
faint flux levels, it becomes ever more difficult to recognise 
sparsely populated stellar systems. Important progress has 
been made recently resulting from newly discovered ultra- 
faint dwarf galaxies, which have nea rly doub led the num- 
ber o f known dwarfs ( Willman et al.li2005 : Belokurov et al.l 
120071 ). Yet the dwarf counts in the Milky Way remains far 
to small compared to the number of subhalos predicted by 
simulations. 

Interestingly, the discovery of compact hydrogen clouds 
potentially located in the Galactic halo has revived the 
possibility of a larger population of gal axy candidates 
that could be traced by neutral hydrogen (JLockmanl l2002l : 
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iRvan- Weber et al.l 200g). This point was made even ear- 
lier by iKlvpin et al.. (,1999. ) who recognised that the much 
broader and physically distinct set of enigmatic HI struc- 
tures commonly referred to as high-velocity clouds (HVCs) 
could represent the missing dark matter subhalo popu- 
lation (IOart]ll970l : iBreemanl Il980l : iBraun fc Burtonll 19991 : 
I Blitz et al .11 19991 ). Although intriguing, to date, a firm HI 
cloud -subhalo link has not been established (JQuilis fc Moord 
l2008l ). 

Prompted by the r ecent release of t he GALFA HI Com- 
pact Cloud Catalogue (jSaul et al.ll2012l ). we revisit the pos- 
sibility that some compact hydrogen clouds can be used as 
a proxy for missing dark matter subhalos. There are inher- 
ent problems when trying to test this possibility, especially 
the lack of distance and mass information for such systems. 
Another major hurdle is the absence of associated stellar 
populations in most compact hydrogen clouds. In fact, it 
is eve n possible tha t many of th e se sys tems never formed 
stars (|Ricottill2009l ). iLewis et al.l (|200O ) proposed exploit- 
ing "pixel gravitational leasing" as a way to map the dark 
matter content in hydrogen clouds. However, the shortage 
of properly aligned nearby background galaxies prevents a 
generalised application to a large sample. 

For this reason, we explore a different approach to de- 
termine their dark matter content. If compact hydrogen 
clouds are dark matter dominated systems, nearby dense 
objects could potentially produce a detectable dark mat- 
ter annihilation s i gnal in gamma rays ( Bergstrom et akl 
1 19991 : iBaltz et al.lbood ). In that vein, iFlix et al.l (|2005l ) 
looked for spatial coincidences between unidentified EGRET 
sources and HVCs. With the vast improvements in gamma- 
ray sensitivity and angular resolution afforded by Fermi 
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Figure 1. Aitoff projection of tlie Fermi all-sky LAT map, showing the locations of the 20 galaxy candidates in our study. Galactic 
coordinates. 



IjAtwood et al.! 120091 ) , we can search for annihilating dark 
matter in GALFA-HI compact hydrogen clouds directly. 
This study should serve not only to investigate a possible 
tracer of gas-bearing dark matter seeds that never formed 
stars, but also to set a general upper bound on the annihila- 
tion cross section {(Jv)^ of a hypothesised weakly interacting 
massive particle (WIMP). Searches for dark matter signals 
using Fermi have been conducted unsuc cessfully in a wide 
variety of astrophysica l systems tBucklev fc Hooper 2010l: 
Ackermann et al.ll20lil : lAndo fc Nagaill2012l : iMirabal et all 



Table 1. Gamma-ray flux upper limits at 95% confidence level. 



2013) • This work simply extends the hunt to a new sample 



The paper is structured as follows. In Section [2] we ex- 
plain the selection of high-latitude galaxy candidates from 
the GALFA-HI Compact Cloud Catalogue, as well as the 
Fermt LAT analysis. In Section|3]we derive gamma-ray flux 
upper limits for Segue 1 and set dark matter annihilation 
constraints for the galaxy candidates. Finally, we briefly 
summarise our interpretation and possible future directions 
in Section m 



GALAXY CANDIDATES AND FERMI LAT 
ANALYSIS 



The GALFA-HI Compact Cloud Catalogue (|Saul et all 
l2012r ) is generated from the Galactic Arecibo L-Band 
Feed Array HI (GALFA-HI) Survey Data Release One 
l|Peek et all 1201 ll ). At completion, GALFA-HI wiU cover 
13,000 deg^ of the sky in the 1420 MHz hyperfine transi- 
tion of hydrogen between Vlsr = ±650 km s""*". Using a 
novel cloud detection algorithm, ISaul et al.l (120121 ) identi- 
fied a total of 1964 compact (< 20') hydrogen clouds in the 
initial 7520 deg^ Data Release One. The catalogue breaks 
down the clouds into a scheme that includes high-velocity 
clouds (HVCs), galaxy candidates, cold low- velocity clouds 
(CLVC), warm low- velocity clouds, and warm positive low- 
velocity clouds in the third Galactic quadrant. For our pur- 
poses, we are only concerned with the 27 possible galaxy can- 
didates that might form the core sample of potentially undis- 
covered dark matter subhalos. In order to guard against pos- 
sible gaseous disk interlopers that may have been pushed 
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prune galaxy candidates at |fe| $; 10°. Finally, we are left 
with a subset of 20 high-latitude galaxy candidates. In Fig. 
[1] we show the distribution of these systems on the sky. 

In order to explore the gamma-ray emission, we use 
the publicly available dataset acquired by the Large Area 
Telescope (LAT) instr ument on board the Fermi Gamma- 
ray Space Telescope (|Atwood et all l2009l 'l. The LAT is a 
pair-conversion gamma-ray detector sensitive to photon en- 
ergies from 20 MeV to 300 GeV. We retrieve all photons 
of 'source' class (evclass=2) within a 10° circular region 
centred at the position of each galaxy candidate. The data 
analysed here were collected between 2008 August 4 and 
2012 November 20 (approximately 4.3 years of data). Good 
time intervals were processed using the available v9r27pl 
Fermi Science Tools with the standard P7S0URCE_V6 instru- 
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Figure 2. Smoothed _Fermi-LAT (1-300 GeV) count map of 
galaxy candidate 143.7+12.9+223. The small white circle marks 
the centre of the Region of Interest (ROI). The map corresponds 
to a 20° circular region. 



ment response function. Throughout, we apply a maximum 
zenith angle cut of 100° . We further filter the data using the 
gtmktime filter expression recommended by the LAT team, 
namely "(DATA_QUAL==1) && (LAT_CONFIG==l) && 
ABS(ROCK_ANGLE)<52" The final analysis for each re- 
gion i ncludes all the poi nt sources listed in the 2FGL cata- 
logue (JNolan et al.ll2012l '). the current Galactic diffuse emis- 
sion model gal_2yearp7v6_v0.f its, and the extragalactic 
isotropic model iso_p7v6source.txt. 

The resulting dataset is analysed with a binned likeli- 
hood method using the gtlike tool in the standard Fermi 
Science Tool^j. For each position, we create a count map 
made up of 30 logarithmically uniform energy bins using 
gtbin. We next construct a binned exposure map with 
gtexpcube2, and a model source and diffuse count map with 
gtsrcMaps. Fig.[2]shows a typical Fermi count map from the 
sample. In the absence of emission, flux upper limits are then 
derived using the implementation of LATAnalysisScriptqj. 
This set of Python libraries unifies the pyLikelihood mod- 
ule included in the standard Fermi Science tools. Upper 
limits are computed with calcUpper assuming a power law 
spectrum of high-energy emission E^^ within a radius of 
10°. We restrict our analysis to photons in the 1-300 GeV 
energy range. The Fermi LAT point spread function (PSF) 
is typically 0.8 (iS/lGeV)""** deg, which in our selected en- 
ergy rang e restricts the photons to less than 1° ar ound each 
location (JGeringer-Sameth fc Koushiappad 120121 ') . Table [T] 
summarises the 95% confidence level upper limits. Since the 
Fermi exposure is rather uniform over the sky the upper flux 
limits are fairly similar, except for a handful of locations with 
higher diffuse background emission or neighbouring bright 
gamma-ray sources. 



3 DARK MATTER CONSTRAINTS 

If compact hydrogen clouds are highly dark matter domi- 
nated objects their gamma-ray emission should be well char- 
acterised by a differential spectrum that can be written as 
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where 
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is the thermally averaged annihilation cross 



section, m^ is the dark matter particle mass, and -77^ is 
the p hoton spectrum of annihilation products (Abd o et al.l 
[2OI0I). The second term, or the so-called astrophysical fac- 
tor J(Ar2), corresponds to the integration of the dark mat- 
ter density squared p^{l,fl) along the line of sight I of the 
compact cloud, over a solid angle AH, so that 
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In principle, we expect compact clouds to be virtually 
free of gamma-ray emission from embedded diffuse and in- 
dividual point sources. As a result, it is reasonable to expect 
that the differential spectrum from dark matter annihilation 
should dominate any gamma-ray signal in the direction of 
observation to these systems. 

Formally, a robust computation of an upper limit on the 
annihilation cross section {crv)~^ requires some knowledge of 
the astrophysical factor J(ASl) of the system under consid- 
eration. Given a lack of direct observational constraints on 
the astrophysical factors of these galaxy candidates Jgc, we 
are only able to estimate (crw)x bounds from our sample 
by tying them to potentially similar systems in the Milky 
Way. Under the Ansatz that they are strongly dark matter 
dominated, the newly discovered ultra- faint dwarf galaxies 
would appe ar to be the closest relatives to compac t hydro- 
g en cl ouds (jStrigari et al.ll200g ). As shown by Sim on et al.l 
(J201H ) , Segue 1 represents the darkest of these systems with 
a very high mass-to- light ratio (~ 3400 Mq/L©) and dark 
matter density 2.5ljg M© pc""^. It also boasts the largest 
astrophysical fac tor for known dwarfs Jseguei = lo'^® 
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Hereafter, we derive upper limits on (crv)^ relative to 
the bounds already imposed for Segue 1 using Fermi mea- 
surements (lAbdo et al.ll2010l: lEssig et al.ll201Gl : IScott etlo] 
bold : iGeringer-Sameth fc Koushiappad I2OI2I ') . Accordingly, 
we repeat the previous Fermi LAT analysis now centred on 
Segue 1 {£,b) = (220.5°, 50.4°). The corresponding gamma- 
ray upper limit for Segue 1 is Fum seguei (1-300 GeV) = 
3.5 xlO~^^ ph cm~^ s~^. In order to turn our flux upper 
limits into a bound on {(Jv)^, we need to assume a spe- 
cific annihilation channel. We adopt the strictest possible 
limit on dark matter annihilation into bottom quarks bh for 
Segue 1 based on 3 years of Fermi data, which translates into 
<^v)y 5. 3xl0~^^ cm^ s~'^ for particle m asses m^ < 40 GeV 
Geringer-Sameth fc Koushiappasll2012r ). Assuming the an- 
nihilation to be purely into bh and an average flux for the 
galaxy candidates {Fum{gc)) — 7.7 x 10"^^ ph cm~^ s~^, 
we can approximate (uv)^ from these systems as 
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for m^ < 40 GeV WIMP masses annihilating to bb, depend- 
ing on the yet unknown Jgc- 



4 DISCUSSION AND FUTURE PROSPECTS 

We have reported gamma-ray flux upper limits for a subset 
of galaxy candidates from the GALFA-HI Compact Cloud 
Catalogue. It is difficult to imagine any of these systems 
beating the astrophysical factors of known dwarf galaxies 
that spans the gamut from 4 x lO^'^ GeV^ cm~^ for Carina to 
1. 3 X 10^^ GeV^ cm ~^ for the ultra-faint dwarf galaxy Segue 
1 (jEssig et al.ll2010l V Without mass and distance constraints, 
it will be difficult to place compact clouds in the context of 
other measurements. It is generally assumed that in order 
to be self gravitating, compact clouds could pack masses 
as low a s a few Mp to greater than ^ 10^ M p for outer 
systems (|Giovanelli et al.ll2010l : ISaul et al.ll2012h . Our naive 
expectation is that Jseguei/Jgc 3> 1 should hold in most 
cases. 

As an illustration, the average compact cloud properties 
indicate a mass of ~ 2 x l O'* M^ and a phy sical size of ~ 100 
pc at distance of 100 kpc IjSaul et al.ll2012l ). Since the astro- 
physical factor is proportional to the density squared, we 
can roughly approximate Jseo uei/Jpc « 10^ for Mseguei ~ 
6 X 10^ Mq (|Simon et al.ll201ll '). To be sure, arguments con- 
sistent with an interpretation that Segue 1 is a tidally dis- 
rupting star cluster contaminated b y the Sagittarius stream 
should be definitively ruled out IjNiederste-Ostholt et al.l 
[2009). For the entire family of ultra- faint galaxies, it is 
also critical to deriv e more reliable estimates of the J-factor 
l|Walker et al.ll201ll ). 

This null result joins the ranks of past dark matter an- 
nihilation searches, which have failed to detect a gamma- 
ray signal in systems suspecte d of high dark matter content 
(jBringmann fc Wenigeill2012l 'l. Unfortunately, as with the 
rest of dark matter pursuits, this is an everything or nothing 
undertaking. Here, we are further hampered by the fact that 
we are seeking a detection with nearly no information about 
distances and masses of the objects involved. Nonetheless, 
Fermi will continue to collect data and stricter gamma-ray 
limits for these systems can be reached. Above E > 100 
GeV, the Cherenkov Telescope Array (CTA) will be crucial 
to escalate the dark matter search to unp recedented bounds 
||CTA Consortiumll201ll : lDoro et al.ll2013l '). There is also suf- 
ficient motivation to explore signatures for ot her reasonab le 
dark matter candidates at other wavelengths (|Feng|[2010|). 

From our measurements, it is still unclear where 
GALFA-HI galaxy candidates fit in the larger dark matter 
subhalo picture. Searches for stellar counterparts associated 
with th ese gas-bearing s ystems are underway and should in- 
tensify (jSaul et al.ll2012h . If this goal is accomplished, mem- 
ber stars could be used to produce a dark matter density 
profile. In an optimistic scenario, compact clouds could val- 
idate a sort of "hiding in plain sight" model whereby dark 
matter subhalos at small scales would be traced directly 
through neutral hydrogen. 

Alternatively, dedicated observational studies could fi- 
nally establish that all compact clo uds are purely bary- 
onic and hence devoid of dark matter IjPlockinger fc Henslen 
[2013). The disagreement between observations and ACDM 
simulations might then have to be invoke more inventive so- 



lutions (JBovlan-Kolchin. Bullock fc KaplinghaJlJOl^ 'l. Re- 
gardless, the potential dark matter fingerprint discussed here 
stands a possible diagnostic of suspected nearby dark matter 
dom inated galactic can didates. With the advent of Skymap- 
per (IKeller et al.ll2007^ a nd the Large Synoptic Survey Tele- 
scope l lvezic et al.l 120081 ) . it might be possible to trace the 
subhalo population directly using faint stars. A cross-match 
between stellar concentrations and compact hydrogen cloud 
positions should be conducted as soon as said surveys are 
completed. 
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